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We address the general problem of hard objects on random lattices, and emphasize the 
crucial role played by the colorability of the lattices to ensure the existence of a crystal- 
lization transition. We first solve explicitly the naive (colorless) random-lattice version of 
the hard-square model and find that the only matter critical point is the non-unitary Lee- 
Yang edge singularity. We then show how to restore the crystallization transition of the 
hard-square model by considering the same model on bicolored random lattices. Solving 
this model exactly, we show moreover that the crystallization transition point lies in the 
universality class of the Ising model coupled to 2D quantum gravity. We finally extend 
our analysis to a new two-particle exclusion model, whose regular lattice version involves 
hard squares of two different sizes. The exact solution of this model on bicolorable random 
lattices displays a phase diagram with two (continuous and discontinuous) crystallization 
transition lines meeting at a higher order critical point, in the universality class of the 
tricritical Ising model coupled to 2D quantum gravity. 
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1. Introduction 



In lattice statistical mechanics, universality classes usually do not depend on the 
lattice over which the model is defined, but only on the symmetries of the interactions. The 
situation however becomes more subtle when the symmetries of the interactions themselves 
strongly depend on the properties of the underlying lattice. A famous example of such 
behavior is provided by the general problem of "hard" objects on two-dimensional lattices. 
In these models, each site of the lattice may be in two states, occupied or empty, but if 
a site is occupied, then necessarily its nearest neighbors must be empty. The models are 
further defined by attaching an activity z per occupied vertex. For regular lattices in two 
dimensions, an exact solution for the thermodynamics of the model exists so far only in the 
case of the triangular lattice (hard hexagon model, solved by Baxter In this case, two 
critical points were found at values z± = (^^^)^, Z- < 0, 2+ > 0, respectively governed 
by the Lee- Yang edge singularity in 2D |^,^ (non-unitary Conformal Field Theory (CFT) 
with central charge c(2, 5) = —22/5) and the critical three-state Potts model (unitary 
CFT with central charge c(5, 6) = 4/5). Despite recent progress the hexagonal 

lattice and square lattice cases remain elusive 0. Numerical evidence however seems to 
indicate that they still have two critical points z- < and z+ > 0, and while the first 
still corresponds to the Lee- Yang edge singularity, the other one displays the exponents 
of the critical Ising model (CFT with c(3, 4) = 1/2). The difference in universality 
class at zj^ for the triangular lattice on one hand and the hexagonal or square lattice 
on the other may be simply understood from the symmetries of the lattices. Indeed, 
this critical point corresponds in all cases to a crystallization transition, where the hard 
particles occupy preferentially a particular sublattice of the lattice at hand. The common 
feature of the square and hexagonal lattice is their vertex-bicolorability (bipartite nature) 
which naturally defines two equivalent mutually excluding sublattices corresponding to two 
possible symmetric crystalline groundstates. The triangular lattice on the other hand is not 
vertex-bicolorable, but vertex-tricolorable instead, which allows to define three equivalent 
mutually excluding sublattices corresponding to three possible symmetric crystalline states. 
These two- or three-fold symmetries give rise naturally to critical points with ^2 (Ising) 
or ^3 (three-state Potts) symmetries. 

The purpose of this note is to study similar hard particle models on random lattices 
such as those used to generate discrete models of 2D quantum gravity. The aim of this 
study is twofold: (1) to check the crucial role played by the vertex-colorability of the 
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underlying lattices in determining the physical behavior of the models and (2) to give a 
"gravitational" proof that the crystallization transition point z+ indeed lies in the critical 
Ising universality class in the case of hard particles on bicolorable lattices. 

More precisely, in this paper, we first solve explicitly the problem of hard particles 
on arbitrary random lattices and show that in the absence of any colorability constraint 
(Euclidean random surfaces), only the Lee- Yang critical point survives at some negative 
value z = Z-. We then solve the same model on so-called Eulerian random surfaces, i.e. 
random lattices for which we impose vertex-bicolorability, and find that the crystallization 
transition point z = is restored in this case, and that it belongs to the universality class 
of the critical Ising model coupled to 2D quantum gravity. 

We next extend our analysis to higher order critical points and show how to recover the 
tricritical Ising universality class (CFT with central charge c(4, 5) = 7/10) by considering 
a two-particle exclusion model on vertex-bicolorable lattices. We again give a gravitational 
proof of this fact by explicitly solving the model on random vertex-bicolorable lattices. 

The paper is organized as follows. In Sect. 2 we study the model of hard particles and 
first recall a few known facts on its regular honeycomb and square lattice version (Sect. 2.1). 
We then solve the model in Sect. 2.2 on arbitrary random planar tetravalent lattices, by 
use of a two-matrix integral. In addition to the activity z per particle, the latter include 
an extra weight g per vertex (occupied or not). We derive explicitly the gravitational 
critical line g = gdz) selecting arbitrarily large lattices, allowing to reach the interesting 
thermodynamic behavior of the model. On this line, we find a unique critical point of the 
matter system at some z = Z-, which we identify as the Lee- Yang edge singularity. In 
Sect. 2.3, we introduce a four-matrix model describing the same problem now on random 
vertex-bicolorable lattices. Wc solve this model exactly in the case of trivalcnt such lattices 
and obtain in particular the new gravitational critical line g = gc{z). Along this line, we 
now find two matter critical points at some values z = Z- < Q and z = Zj^ > The first 
one is still identified as the Lee- Yang edge singularity, while the new one is identified as 
a crystallization transition point in the universality class of the critical Ising model. Sect. 
3 is devoted to the study of a more sophisticated two-particle exclusion model in which 
we allow sites to be occupied by single particles (with activity zi per particle) or pairs 
of particles (with activity Z2 per pair), while the exclusion rule imposes that each edge 
of the lattice is shared by a total of at most two particles. In Sect. 3.1, we present its 
expected phase diagram on a regular square or honeycomb lattice obtained by applying 
the same ideas as above, i.e. relating the different phases of the model to the compatibility 
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between the exclusion rules on one hand and the vertex-bicolorability of the lattice on the 
other. These ideas are tested in Sect. 3.2 where we solve the model on random trivalent 
vertex-bicolorable lattices, by means of a six-matrix integral. The gravitational critical 
surface g = gc{zi,Z2) is explicitly shown to contain the expected matter phase diagram, 
with in particular a critical Ising transition line meeting a first order line at a tricritical 
Ising transition point for some positive values of the activities (21,^2) = {z^^\z2^). It 
also displays a line of Lee- Yang edge singularities terminating at a higher order critical 

point described by a non-unitary CFT with central charge c(2, 7) = —68/7 at some point 

it') it') 

zi = z\ > and Z2 = Z2 < 0. A few concluding remarks are gathered in Sect. 4, while 
additional technical derivations or more involved cases are left to appendices A-E. 



2. Nearest neighbor exclusion models on regular and random lattices 

In this section, we study generic models of nearest neighbor exclusion in which hard 
particles live on the vertices of either regular on random lattices. The exclusion rule 
simply states that when a vertex is occupied by a particle, all its nearest neighbors must 
be vacant. We first briefly recall known facts about the so-called "hard square" model 
of hard particles on the regular square lattice [|^-|^ and the corresponding model of hard 
particles on the regular honeycomb lattice ||^ . As explained below, these models share the 
same qualitative phase diagram, where the structure of the ordered phase strongly relies 
on the vertex-bicolorability (i.e. the bipartite nature) of the underlying lattice. We then 
present an exact solution of the same models on random lattices. We first consider the 
case of arbitrary random graphs standardly used in the context of discretized 2D quantum 
gravity P,p!0[]. We note the disappearance of the ordered phase which can be traced back 
to the generic lack of vertex bicolorability of the graphs. We therefore study the case of 
vertex-bicolorable (or so-called Eulerian) random graphs, for which we derive the phase 
diagram in the planar limit. As expected, the ordered phase is reinstated in this case, and 
we find a continuous transition point in the universality class of the critical Ising model 
coupled to 2D quantum gravity. 
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Fig. 1: Sample configurations of tlie liard particle model on the honeycomb 
(a) and square (b) lattices. The exclusion constraint between particles (rep- 
resented as black dots) may be translated into a non-overlapping constraint 
for hard tiles, either hexagons (a) or squares (b). 

2.1. Nearest neighbor exclusion on the square and honeycomb lattices 

Let us consider nearest neighbor particle exclusion models on respectively the regular 
trivalent honeycomb lattice and regular tetravalent square lattice. Each particle comes 
with an activity z and excludes its nearest neighbors. A simple pictorial representation 
of the exclusion rule is to replace particles with non-overlapping hexagons, resp. squares, 
with centers on the vertices of the lattice, as depicted in Fig. 1. The phase diagrams of 
both models have been derived using numerical (corner) transfer matrix methods in 
They both display two critical points. 

The first "non-physical" critical point occurs at a negative value Z- of the activity 
and corresponds to the so-called Lee- Yang edge singularity |jTl|] (non-unitary CFT with 
central charge c(2,5) = —22/5), with the estimates z^'* = —0.15 . . . (extracted from 0) 
for the (trivalent) honeycomb lattice and z^_^^ = —0.122 . . . for the (tetravalent) square 
lattice. The values Z- of z at these critical points can be obtained from the singularity 
of the thermodynamic free energy F{z) expressed as an alternating series in powers of z. 
For z Z-, the singular part of the free energy behaves as F{z)\sing ~ (^ — Z-)'^~" where 
a = 7/6 is the thermal exponent predicted by the CFT. This apparently non- physical 
critical point with negative activity can be reinterpreted as a positive activity critical 
point at = —Z- for heaps of hexagons, resp. squares with activity t = —z per object in 
2-1-1 dimensions, with free energy $(t) = —F[z = —t) [ p!2| . 

The second "physical" critical point occurs at a positive value z+ and corresponds to 
a continuous transition between a low activity disordered fluid phase and a high activity 
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ordered crystalline phase where the hard particles condense preferably on one of the two 
sublattices of the lattice. Estimates for the critical points are z\_ = 7.92 . . . [|| and 
z^'^ = 3.7962... 0. More precisely, the honeycomb and square lattice are bipartite 
lattices whose vertices can be naturally bicolored, say black and white, in such a way that 
a vertex of one color has only neighbors of the other. The corresponding order parameter 
M = pb — Pw-, which measures the difference of density of particles between the black 
and white sublattice, is zero in the fluid phase and non-zero in the crystalline one. Based 
on numerical evidence, it is commonly accepted that this transition lies in the universality 
class of the critical Ising model (CFT with c(3, 4) = 1/2), with a thermal exponent a = 
and a magnetic exponent /3 = 1/8 corresponding to the singularity of the order parameter 
M ~ (2 — Zj^Y z ^ 

Note that the two above critical points are essentially different in nature. The Lee- 
Yang critical point is much more universal and does not rely on any particular geometrical 
feature of the lattice. It is also observed for hard particles on the triangular lattice (the 
so-called "hard-hexagon" model). On random lattices, it was first observed for the Hard 



Dimer model in [|T3[ and, as we shall see, it will show in all the hard particle models 
studied throughout this paper. On the contrary, the universality of the crystallization point 
strongly relies on the bicolorability of the underlying lattice. In the case of the triangular 
lattice for instance, one finds instead a crystallization point governed by the critical three- 
state Potts model (CFT with c(5,6) = 4/5), directly related to the tricolorability of the 
lattice. In the case of random graphs, we shall also find that the bicolorability is crucial 
to recover the Ising crystallization point. 

2.2. Nearest neighbor exclusion model on a random lattice 

We now turn to the study of nearest neighbor exclusion models on random graphs. 
In this section we will concentrate on the case of hard particles living at the vertices 
of arbitrary random tetravalent graphs. We use the standard matrix integral method to 
generate all configurations of the exclusion model on all possible tetravalent fatgraphs, 
including a weight g per tetravalent vertex (empty or occupied), z per particle, and the 
usual weight A^2-2/i graphs of genus h. We will mainly be interested in the planar limit 
N —>■ 00 that selects only graphs with the topology of the sphere h = 0. More precisely, 
the generating function reads 



(2-1) 

V{AB) = ~-A'+AB-g—-gz— 



where A, B are Hermitian matrices with size N x N, and the measure is normahzed so that 
Z^\o, 0) = 1. The Feynman diagrammatic expansion of (|2.1| ) is readily seen to generate 
tetravalent graphs with two types of vertices, occupied with weight gz) and empty 

with weight ^r), and with the propagators given by the inverse of the quadratic part 
of V, namely: 

l',»a = i(-4,B)Q(^) and ( <^^> ]) P-^) 

hence {BijBki) = dudjk/N and {AijBki) = dudjk/N (for a general review on matrix models 
see P,pil[] and references therein). The vanishing of the {AA) propagator clearly enforces 
the exclusion rule. 

To compute (^T|) 5 ripply the standard technique of bi-orthogonal polynomials |§] . 
We introduce the monic bi-orthogonal polynomials pn, Qm with respect to the scalar prod- 
uct (f^g) = J dxdye~^^^^'y^ f{x)g{y), satisfying 

where the "norms" are fully determined by the requirement that the polynomial Pn be 
monic of degree n. Introducing hn ^ = h^ig = 0, z = 0), and after reduction of ( |2.1| ) to an 
eigenvalue integral, we may rewrite 

ZNi9,z)=ll^ (2.4) 

n=0 '^■n 

which reduces the computation of the partition function to that of the h^^s. 

We introduce the operators of multiplication by eigenvalues Qi and expressed 
on the orthogonal polynomials as QiPn{x) = xpn{x) and Q2<lm{y) = yQmiy)^ and the 
operators of derivation with respect to eigenvalues Pi and P2 expressed as PiPn{x) = p'^ix), 
P2(lm{y) = Qmiy)- Integrating by parts, we get the system of equations 

d 

[PiPn.qm) = N{—V{x,y)pn{x),q^{y)) = N{pn,Q2qm) - N{{Qi + gzQi)pn,qm) 

Q (2.5) 

{Pn,P2qm) = N{ — V{x,y)pn{x),qm{y)) = N{QiPn,qm) - Ng{pn,Qlqm) 
dy 

Note that QiPn is a linear combination of the pj^s for < j < n -|- 1, and similarly for 
combination of the g^'s for < z < m -|- 1, while PiPn is a linear combination of 
the Pj 's for < j < n — 1 and similarly P2Qm, combination of the g^'s for < z < m — 1. 



6 



Prom orthogonality, we see that for m > n + 3 in the first hne of ( |2.5|) and for m < n — 3 
in the second, we get respectively that (pn, Q2?m) = and {QiPmQm) = 0. Hence the 
linear combinations reduce respectively to finite ranges of indices n — 3<j<n + l and 
m — 3<i<m + l. Moreover, from the parity of the potential V{x, y) = V{—x, —y), we 
deduce the parity of the polynomials: Pn{x) = {—l)'^pn{—x) and qm{y) = (— l)™'?m(— ?/)■ 
So finally the action of Qi and Q2 takes the form 

QlPnix) = Pn+l{x) + rnPn-l{x) + SnPn-six) 

(2.6) 

^mQm— 1 {y) + {y) 

The fact that Qi and Q2 have a finite range is generic of multimatrix models with polyno- 
mial interactions. Eqns. (|2.5| ) and ( |2.6|) can be expressed in an operatorial way. Introducing 
the adjoint operators Q\ and Q2 with respect to the above scalar product, eqn. ( |2.5| ) takes 
the form: 

Ql-Qi- gzQt 



^ (2.7) 

-^ = Q{-9Ql 

Let us introduce the shift operators cr, r acting respectively on the p's and g's as apn = 
Pn+i, Tqn = Qn+i and their adjoints cr^,r^, such that 

a^ = T~^v T^ = a~^v (2.8) 

where v = is the diagonal operator acting as vpn = v^Pn and vq^ = VnQm with 



hn-l 



(2.9) 



Analogously we define the diagonal operators i^, r, s, r, 5 acting on p^ and respectively as 
the multiplication by n, Sn,fn,Sn- In terms of these operators, the P and Q operators 
read finally 

Qi=a + a~^r + a~^s 

Q2=T + T~^f + T~^S 

q\ = r~^v + rv~^T + s{v~^Tf 

(2.10) 

Q\ = a~^v + rv~^a + s{v~^af 

= (j-V + 0((7-2) 
P2 = r-V + 0(r-3) 
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To obtain a system of recursion relations involving the sequence = hn/hn-i, let us now 
write order by order in r and a the two relations (|2Tj 



0(a"^) : a"^-^ = a ^{v -r) - gz{ra V + a + (a ^rfa + a ^s + a '^sa + a ^sa^) 
0{a) : = fv~^a — a — gz{ar + ra + a~^ra'^) 
0{a^): = s{v-^af - gza^ 

(2.11) 

and 

0(r) : = rv~^T — g{Tf + rr + T~^fT'^) 
0{t^): = siv-^rf - gr^ 

(2.12) 

It is not difficult to check that the first lines of (2.11) and (2.12) are equivalent modulo 



the other equations, so that we are left with five equations for the five unknown sequences 

^715 ^n- 

The planar (genus zero) limit of these equations amounts to taking n, N ^ oo with 
X = n/N fixed, in which case all sequences converge to functions of x. More precisely, we de- 
fine the rescaled limits V{x), i?(x), S'(x), -R(x), S{x) of respectively gvn, grn, g^Sn, gfm g'^Sn- 
These functions are determined by rewriting ( |2.11| ) and ( |2.12|) in this limit which amounts 
to treating all operators as scalar s (in particular a = r = 1), with the result 

gx = V-R- 3z{R^ + S) 
^ R-V(l + 3zR) 

= S-zV^ (2.13) 
= R~ 3VR 

o = s-v^ 

After substitutions, this reduces to 

gx = ^{V) ^ V{1 - 3zV^) - (2.14) 

For fixed g and z this equation defines upon inversion the function V{x) encoding the 
asymptotic properties of the sequence f^, hence those of the h^s. More precisely, using 



eqn. (|2.4|) , the thermodynamic free energy in the planar hmit reads 

/(^)(r7, z) = - lim Log(4^)(^7, z))/N^ 

N^oo 



ax(l — xjLog 



\ 

where the normahzation gx in the logarithm ensures the correct normalization of the 
partition function, namely fQ^\g = 0, ^ = 0) =0. For g sufficiently small, eqn. ( |2.14| ) 
determines a unique solution V{x) which is monotonous and such that V{x) ~ gx for 
small X. To compute fQ^\g,z), we must substitute this solution into ( |2.15|) and perform 
the integration. It is natural to perform the change of variables x ^ V, which yields 

where Vg^z is the value of F at x = 1, satisfying Lp{Vg^z) = 9 for fixed g and z. The 
singularities of the planar free energy are due to those of Vg^2 as a function of g and z. 
For fixed z, the first singularity of fQ^\g: z) is attained at a critical value g = gc{z), where 
the value Vc{z) = Vg^(^z),z is such that {p'iVc^z)) = 0. In the vicinity of this point, ( |2.14| ) 
reduces to gc{z) —g ~ iVc^z) — Vg^zY, which in turn yields a generic square root singularity 
for Vg^z in (|2.15| ). To get the corresponding singularity of the free energy, we note that, 
taking successive derivatives of ( p.l6|) with respect to g and using ^{Vg^z) = 9-, we simply 
get the general formula 

f.o^Tf'^'M-^-^'^ (2.17) 

dg^ dg Vg^z dg 

The square root singularity of Vg^z immediately translates into a string susceptibility ex- 
ponent 7, defined by /o^''(5', -2)|sing ~ {9c{z) — gY~"' , with value 7 = —1/2. The critical 
value gc{z) also corresponds to the maximum value of g for which the series expansion of 
the free energy converges. 

Writing gc{z) = ifi{yc{z)) and = (p' {Vc^z)) yields parametric equations for the critical 
line gc{z), in terms, say, of the parameter u = 3zVc{z)'^: 

_12u{l + 3u)^ ^ ^ _ (1- 3^)^(1 + 10^-15^^) 



(l-3w)8 ' ' 12(1 + 3m) 



This line is plotted in Fig. 2. The solid curve corresponds to the true critical values of 
9 = 9c{z), corresponding to a first maximum of ^{V), while the dashed line corresponds 
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Fig. 2: Critical line gdz) for the hard particle model on arbitrary random 
lattices, as obtained by setting ip'{V) = 0. The solid line represents the 
true critical values of g = gdz) corresponding to the maximum of (f first 
attained by the change of variable x ^ V, while the dashed line corresponds 
to a further minimum as displayed on the top right plot of the figure. These 
extrema merge at the Lee- Yang critical point Z- where f"{V) = (top center 
plot) and disappear from the real plane (by becoming complex) for z < Z- 
(top left plot). 

to a further minimum of f{V) which is never attained by the above change of variables. 
The solid line stops at a finite negative value of z below which gdz) becomes complex. 
We note that for z = {u = 0) we recover the critical value gc = 1/12 for pure arbitrary 
tetravalent graphs WM- The solid and dashed lines merge into a cusp at a higher order 



critical point satisfying in addition f"{V) = 0, as it corresponds to the coalescence of the 
maximum and the minimum of ip (see Fig. 2). This point corresponds to the values 

25 



8192 



(11^5 + 25) = -0.151... 



g- = gdz-) = ^(13v^ " 29) = 0.0979... (2-19) 
V- = Vdz-) = 32/75(7v^- 15) = 0.278... 

At this point, the above scaling argument now becomes g- — g ^ {V- — V)^, hence 
translates directly into a string susceptibility exponent 7 = — 1/3. As the critical activity 
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Z- is negative, we identify this higher critical point with the unique non-unitary theory 
with 7 = —1/3, namely the Lee- Yang edge singularity with c(2,5) = —22/5, coupled to 
2D quantum gravity (see Appendix A). This identification may be corroborated in two 
different ways. On one hand, one may compute the thermal exponent a. which measures 
the singularity of the free energy at the critical point as a function of and compare it 
to the predicted value a = 1/2 for the Lee- Yang edge singularity coupled to 2D quantum 
gravity, as re-derived in Appendix A below. Writing the thermodynamic free energy per 
site as f^'^^z) = — lim^^oo \^og Za{z) where Za{z) is the partition function for planar 
graphs of fixed number of vertices (area) A, we read a from the singular part /'■^•'(^)lsing ~ 
{z — z-Y~'^ . Using f^'^\z) = Loggdz), a may be obtained by computing the singular 
part of gc{z). Expanding eqn. ( |2.18|) in powers of w — m_, we find that gdz) — g- = 



a{u — U-)'^ + b{u — U-)^ + ... and z — Z- = a' {u — U-)'^ + b' {u — u_)^ + . . . with U- = 3z-V^, 
and which, upon inversion, yields gdz) — g- = a" {z — Z-) + b" {z — Z-)^ + ... with a 
non- vanishing value of b" . Hence we get 2 — a = 3/2as expected. On the other hand, one 
may also derive the so-called double scaling limit of the model at the critical point and 
write a differential equation for the renormalized string susceptibility, easily identified with 
that of the Lee- Yang edge singularity coupled to 2D quantum gravity. This derivation is 
presented in Appendix B below. 

Note that for z < Z-, the gravitational critical value of g becomes complex with 
gdz) = pe*^, but it still governs the large area behavior of the above partition function 
Za{z) which now oscillates typically as 

ZAiz) p-^cos{Ae) (2.20) 

which allows to identify the thermodynamic free energy as Log p. 

In conclusion, when comparing with the regular lattice results, we see that the naive 
gravitational version of the exclusion model fails to reproduce the crystallization transition, 
and leaves us only with the "non-physical" Lee- Yang edge singularity. As shown in Ap- 
pendix C, the case of trivalent graphs instead of tetravalent is exactly solvable as well and 
displays the same structure. The absence of a crystalline ordered phase should not come 
as a surprise, as the partition function involves a sum over graphs that are not generically 
bicolorable, hence do not allow for a canonical crystalline order, where half of the vertices 
are preferentially occupied. To emphasize the role played by bicolorability for exclusion 
models, we note that for large z {u ^ 1/3 in eqn. ( p.l8| )) the quantity y/z g^z) tends to 
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2/9, which is precisely the critical value of (7 in a model of pure bicolorable tetravalent 
graphsS. This clearly shows that in this limit the selected configurations are half-occupied 
vertex-bicolorable graphs. 

To recover a crystallization transition at finite 2;, we shall consider in the next section 
the coupling of exclusion models to Eulerian gravity. 



2.3. Nearest neighbor exclusion models on vertex-bicolorable random lattices 

In this section, we consider a restricted gravitational version of the nearest neighbor 
exclusion model, in which we explicitly sum only over the so-called random Eulerian graphs, 
simply defined as vertex-bicolorable (or bipartite) graphs. It turns out that the case of 
trivalent Eulerian graphs is technically simpler than that of tetravalent ones, yet it displays 
the same qualitative physical behavior. Therefore we will now concentrate on the trivalent 
case and leave the tetravalent one to Appendix E. 

The configurations of the nearest neighbor exclusion model on Eulerian trivalent 
graphs are again generated by a matrix model replacing (|2.1|) . We now need a total 
of four matrices, as the vertices must be bicolored and empty or occupied. More precisely, 
we use a matrix Ai (resp. A4) for empty black (resp. white) vertices and a matrix A2 
(resp. A3) for occupied white (resp. black) vertices. The resulting matrix model reads 



Z^j^\g,z) = J dAidA2dAsdA4e-''^'^^'^^'^^'^'^^^ 



.3 (2-21) 

ViA^, A2, A3, A4) = A^A2 - A2A3 + A3A4 - gC-^ + ^-^) - gzC-^ + ^ 



.A? A3 Al Ai 



The quadratic form in V{Ai, A2, A3, A4) has been engineered so as to reproduce the correct 
propagators, namely that only black and white vertices are connected in the Feynman 
diagrams ((A^Aj) = if i and j have the same parity) and that two occupied neighboring 
vertices exclude one-another ((A2A3) = 0). 

Due to the chain-like interaction between the matrices, this model turns out to be 
solvable by means of bi-orthogonal polynomials. In addition the symmetry A^ ^ A^^i 
implies that the two sets of polynomials are identical. We therefore introduce the monic 



This critical value was computed in |15] in the context of a particular 0(n=l) gravitational 
model. 
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polynomials Pn, orthogonal with respect to the appropriate symmetric scalar product, 
namely 

{Pn,Pm) = j dXidX2dX3dX4e~^^^''^'''^^''^^''^^PniXl)Pmix4) = hr,5r,,m (2.22) 

As before, we introduce operators Qi of multiplication by Xi, i = 1, 2, 3, 4 but this time all 
acting on Pn{xi), but the symmetry Ai ^ A^^i immediately implies that Qs = Q\ 
Q4 = Q\- We also introduce the operator Pi acting on as d/dxi. These operators 

satisfy the system 

N (2.23) 

Qi = Q\+gzQl 

obtained by integrating by parts. 

To write explicitly the action of Qi and Q2 on the Pn's, let us first notice that the 
potential V satisfies the symmetry relation 

u;^X2, wxs, Ci;^X4) = y(xi, X2, X3, X4), w = e^*'^/^ (2.24) 

This translates into the symmetry relation for the monic orthogonal polynomials 

cc;2>„(wx) =pn(x) (2.25) 

Now from (12^231) it is easy to show that Qi and Q2 have finite range, and more precisely, 
thanks to the symmetry relation ( p.25| ) 

(2.26) 

where a is the shift operator acting on the p's as apn = Pn+i and the operators s^^^ 
are diagonal. Introducing as before the diagonal operator v with entries Vn = /in/^n-i, 
we have = a~^v, and therefore 

gt = s(0)(^-i^)2 + ^(D^-V + s(2)(t;-V)4 (2.27) 
Expanding the relations (|2.23|) order by order in cr, we finally arrive at 

= 9 

V = s^^) + g'^zv{as^^^a~'^ + a'^s^^^a) (2.28) 

Ji = s(i)-(7(r(i) + a-VWa) 
N ^ ' 
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where v is the diagonal operator with entries n. 

In the planar limit, we express the system (|2.28|) in terms of the rescaled limiting 
functions S{x) = lim (/^zs*-^-*, V{x) = Xmig^zv^ where x = n/N while n, N ^ oo, and also 



use the fact that a 1. The third line of (|2.28|) allows to solve for S = V/{1 + 2V), so 
that we finally get 



V 



1 + 2^)2 



21/2(1 _ 21/2) 



Writing 



(2.29) 



(2.30) 



(1 + 2^)3 

we get two candidates for the critical line, namely V = 1/2 or z = 41/(1 + 21/)^. 

These lead respectively to the two possible critical curves g^{z) = Lp{V)/z^ parametrized 

by 



1 1 



9t{z) = ^--rT^ '2>0' V = l/2 



8z 4^2 

2,^, 1 + 8^ + 10^ 2 ^^^^^ ^^^^^ _ ^^^^^ ^^^^^ 

i + 2vy 



with z{V) = 41/(1 + 2Vy 



(2.31) 



where the condition 2; > in the first line simply follows from the positivity requirement 
for the norm ratios f^, implying that V and z have the same sign. 




V 



Fig. 4: Critical lines in the iV^z) plane as obtained by setting ^' iV) = 0. 
The correct line corresponds to the lowest value of |1/| and is represented by 
a solid line. The Lee- Yang critical point Z- is characterized by dz/dV = 
and corresponds to the merging and annihilation of two extrema. The Ising 
critical point corresponds to the crossing of two determinations of V{z). 
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Fig. 3: Critical line gdz) for the hard particle model on vertex bicolorable 
random lattices, as obtained by setting f'{V) = 0. The solid line represents 
the true critical values of g = gdz) corresponding to the maximum of f first 
attained by the change of variable x V, while the dashed line corresponds 
to a further minimum as displayed on the top plots. The curve terminates 
at a Lee- Yang critical point similar to that of Fig. 2. The maximum and 
minimum swap determinations at the critical Ising point z.^. where f'^V) = 
(top center plot). 



These two lines are represented in Fig. 3. The choice of the correct determination 
m ( |2:31| ) is best seen by plotting the critical lines in a {V, z) diagram as in Fig. 4. In 
this picture, the correct singularity of the free energy at fixed z is always given by the 
lowest critical value of \V\ as it is the one attained by the above change of variables. 
The true critical points of the free energy are represented by solid pieces of curves in 
Fig. 3, and correspond to maxima of ^{V), while the dashed parts correspond to further 
minima of ^{V) never attained by the change of variables as before. Along these lines 
the string susceptibility exponent is 7 = —1/2. We note that for 2; = 0, we recover the 
critical value (7^ = 1/8 for pure bicolored trivalent graphs |jl^,0. Also for 2 ^ 00 we get 
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g'^{z) = gf{z) ~ 1/(8^) as expected for half-occupied trivalent bicolored graphs. 

The higher order critical points correspond to the cuspidal singularity at Z-, where 
the maximum and minimum of f{V) merge and annihilate each other, and to the crossing 
between the two curves at z+, where the value of the maximum of f hops from g2iz) to 
gi{z) while its minimum hops from gi{z) to g2{z). The position of these points is obtained 
by writing the extra condition f"{V) = 0, with 



32-z 



if V = l/2 



'^"(^) - ^ 4(21/ -1X1 + lOF) if ;z = 41/(1 + 2y)4 ^^'^^^ 



We get the two critical values 



1 15 

(1) V+ = -, z+ = 32, 9l = ^ 

(2) V- = - — , z- = gt - 



(2.33) 



10' 55' 220 

The critical point Z- is similar to that of Sect. 2.2 and corresponds to the Lee- Yang edge 
singularity. The crucial outcome of our calculation is the emergence of a crystalline phase at 
finite values of z, with a crystallization transition at zj^. This critical point turns out to be 
in the universality class of the critical Ising model on random graphs. The simplest reason 
for these identifications is that both above critical points have string susceptibility exponent 
7 = —1/3 by construction, and that only two CFT's are candidates to describe thisi, the 
Lee- Yang edge singularity with c(2, 5) = —22/5 which is non-unitary as expected for a 
negative value z = Z-, and the critical Ising model with c(3,4) = 1/2 which is a unitary 
theory, as expected for a positive critical activity To further confirm the identification 
of the critical Ising universality class, one can compute the thermal exponent a for the 
crystallization transition. It is easy to see from the above formulas that the transition from 
the curve gi{z) to g2{z) is continuous at zj^, with continuous first and second derivative, 
and with a discontinuity of the third one. This gives a thermal exponent a = —1, as 
expected for the critical Ising model coupled 2D quantum gravity (see Appendix A). We 
have also derived the differential equations for the string susceptibility for both cases in 
the corresponding double scaling limits, and identified them with the known results for the 
Lee- Yang and Ising critical models coupled to 2D quantum gravity. The details of these 
calculations are given in Appendix D below. 



5 Indeed for a CFT with central charge c{p,q) = I — G{p — q)^ / (pq) , we have 7 = —2/{p + q—l) 
hence here p + q = 7, and (p, q) = (2, 5), (3, 4) only. 
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3. Two-particle exclusion models on regular and random lattices 

In this section, we extend our analysis of exclusion models on bicolorable lattices to 
incorporate the physics of higher order critical points. More precisely, in the previous 
section we have reproduced the universality class of the critical Ising model within the 
context of exclusion models on bicolorable random graphs. The crucial feature leading to 
the Ising symmetry is the existence of two degenerate symmetric crystalline groundstates 
playing the role of the two (up and down) ferromagnetic groundstates. 

We now wish to construct in the language of particle exclusion a model reproducing 
the physics of the tricritical Ising model. In the framework of spin systems, the latter is 



found for instance in the phase diagram of a dilute Ising model, with spins a = 0, ±1 118 
where the spins a = play the role of annealed non-magnetic vacancies. The existence 
of a tricritical point is associated to that of three non-symmetric groundstates of constant 
spin: the two groundstates a = +1 and a = —1 play symmetric roles, but a = is on 
a different footing. The Ising second order phase transition in the absence of vacancies 
extends into a line of second order transition points for low enough activity per vacancy. 
On the other hand, at low enough temperature, a first order transition line separates 
a ferromagnetic phase at low activity per vacancy from a paramagnetic phase at large 
activity per vacancy. Both lines merge at a tricritical point whose behavior defines the 
tricritical Ising universality class (CFT with c(4, 5) = 7/10). (For a general review on 
tricritical points in the context of spin systems or lattice gases see fT^.) 

In the next section, we show how to realize a similar behavior within the framework of 
hard objects, by use of a two-particle nearest neighbor exclusion model. We first describe 
its expected phase diagram on a regular bicolorable lattice. Here again, the bicolorability 
of the underlying lattice is crucial for the existence of an ordered crystalline phase. We 
then derive an exact solution of the same model on bicolorable random lattices and recover 
the expected phase diagram, with in particular a tricritical point coupled to 2D quantum 
gravity. 

3.1. Two-particle nearest neighbor exclusion models on the square and honeycomb lattices 

In its simplest formulation, the two-particle exclusion model we wish to study is 
defined by putting particles on the vertices of a regular bicolorable lattice (typically square 
or honeycomb), with the exclusion rule that a total of at most two particles may occupy 
the two vertices adjacent to any edge. In particular, we allow for two particles to occupy 
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(a) (b) 

Fig. 5: Sample configurations of tfie two-particle exclusion model on the 
honeycomb (a) and square (b) lattices. The black dots represent singly occu- 
pied vertices while the circled black dots represent doubly occupied vertices. 
The two-particle exclusion constraint that an edge be shared by at most two 
particles may be translated into a non-overlapping constraint for hard tiles, 
either hexagons for double occupancy and triangles for single occupancy (a) 
or big tilted squares for double occupancy and small squares for single occu- 
pancy (b). 

the same vertex, in which case its nearest neighbors must be vacant, while singly occupied 
vertices may be nearest neighbors. We assign a weight zi for singly occupied vertices and 
a weight Z2 for doubly occupied ones, with typically u = Z2/z1 measuring the attractive 
{u > 1) or repulsive (u < 1) interaction between particles on the same vertex. As in 
the one-particle models of Sect. 2, we may represent pictorially the exclusion constraint in 
terms of non-overlapping hard objects on the corresponding lattices. As shown in Fig. 5, for 
the square lattice, doubly occupied sites may be represented as usual hard squares, which 
exclude their four neighbors, while singly occupied sites are represented by a new type of 
squares which are twice as small and may thus occupy two neighboring sites. Similarly, 
for the honeycomb lattice, we represent doubly and singly occupied vertices respectively 
by hard hexagons and triangles that are twice as small. 

In both cases, we have three crystalline groundstates corresponding to a maximal 
covering of the lattice: two of them are symmetric and use only the larger tiles which 
occupy one of the two sublattices of the bipartite lattice, the third one is obtained by tiling 
the lattice with the smaller objects. The phase diagram of the model is best represented 
in the variables {l/z2,u~^ = zf/z2) (see Fig. 6). On the axis = (i.e. zi = 0), 
we recover the one-particle exclusion models of Sect. 2 with z = Z2, and in particular a 



18 




Fig. 6: Expected phase diagram in the (1/2:2, = z1/z2) plane for positive 

activities. The ordered phase M ^ corresponds to a crystalhzation of the 
doubly occupied vertices on one of the two mutually excluding sub-lattices. 
This phase is separated from the fluid phase M = by a second order critical 
line for small enough and by a first order transition line for small enough 
1/^2 • Both lines are expected to meet at a tricritical point. The three natural 
groundstates in the problem are recovered along the axis I/22 = 0: for < 
1 the two maximally filled doubly occupied configurations dominate (tilings 
with the bigger tiles) while at -u"^ > 1, the full occupation by single particles 
(tiling with smaller tiles) dominates. The latter groundstate degenerates into 
a disordered fiuid phase as soon as 1/^2 > 0. 

crystallization point at Z2 = z+. For fixed small enough u~^, we expect a similar second 
order transition at some Z2 = z+{u). The critical curve Z2 = z+{u) separates the liquid 
phase from the crystal phase, the latter being characterized by the non-vanishing of the 
order parameter M = — expressing the difference of densities of doubly occupied 
sites of either color. On the other hand, on the axis {I/Z2) = 0, i.e. 2:2 — > 00 and zi — > 00, 
with = z\l Z2 fixed, we have a competition between the three groundstates of maximal 
occupation. The two symmetric groundstates made of larger tiles have free energy per site 
(Log 2:2)72, while the other groundstate made of smaller tiles has free energy Log 21. This 
leads to a first order transition at tt = 1, with same order parameter M = ±1 for w > 1 
and M = for tt < 1. We expect this transition point to extend into a first order transition 
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curve for small enough (1/2:2). By analogy with the tricritical Ising model phase diagram, 
we expect the two curves to meet at a tricritical point with c(4, 5) = 7/10. 

These models do not seem to be simply solvable by integrable techniques, and the 
above phase diagram is somewhat conjectural. However, in the following section we shall 
present an exact solution of the random bicolorable lattice version of the two-particle 
exclusion model, and show that it displays precisely the physical picture described above. 

3.2. Two-particle exclusion model on vertex-bicolorable random lattices 




Fig. 7: The non-vanishing propagators corresponding to the six-matrix 
model generating the configurations of the two-particle exclusion model on 
vertex bicolorable random lattices. Identifying respectively Ai, As, A5 with 
empty, singly occupied and doubly occupied black vertices on one hand and 
Aq, A4, A2 with empty, singly occupied and doubly occupied white vertices, 
the selected propagators enforce both the bicoloring constraint (black ver- 
tices are connected to white vertices only and conversally) and the exclusion 
constraint (at most two particles may share the same edge). 

We start with the matrix model 

ZN{g,Zi,Z2) 

V{Ai,A2,A3,A4,A5,Ae) 

where as for Sect. 2.3, the N x N Hermitian matrices Ai with odd (resp. even) index 
correspond to black (resp. white) vertices, which can be empty (^1,^6)5 singly occupied 
(^3,^4) or doubly occupied (^5,^2)- It is easy to check that the inverse of the quadratic 



-NTrViAi,A2,A3,A4,A5,A6) 



i=l 



A1A2 - + - + A^Aq 

A2. 43 /l3 43 43 43 

^4 + t)-^-(T^t)-^-(T^T) 



(3.1) 
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form in V generates the expected non-vanishing propagators (A1A2), {A1A4), {AiAq), 
(AsAg), {AsAq), (743^44) (see Fig. 7). Remarkably enough, the matrix interaction in V is 
simply chain-like, allowing for a solution using bi-orthogonal polynomials. 

As before, the symmetry Ai ^ A-j-i oi V [i = 1,2,3) ensures that the left and 
right polynomials are identical, hence we define the monic orthogonal polynomials Pn with 
respect to the appropriate symmetric scalar product, namely 

6 

^ i=l 

Introducing again the operators of multiplication by eigenvalues Qi, Q3, Q4 = QI, 
Q5 = QI^ Qg = Qi and the operator Pi of derivation with respect to eigenvalues of Ai 
(all regarded as acting on pn{xi)), we find the master equations 

^ = Q2-9Ql 

= Qi-Qs- 9Z2QI (3.3) 

= -Q2 + - gziQl 

which determine the /I's completely. 

As a direct consequence of ( p.3|) , the Q^'s have a finite expansion in powers of the shift 
operator a {apn = Pn+i)- The symmetry relation 

V{uJXi, uj'^X2, uJXs, uP'X/i, wxs, uP'Xf^ = V{xi, X2, X3, X4, X5, xq), uj — e^*'^/"^ (3.4) 

analogous to the case of hard objects on bicolorable trivalent graphs, ensures that the 



relation (|2.25|) still holds for the orthogonal polynomials at hand. Then, using relations 
( p73|) and cr^ = a~^v, where v is still defined as the diagonal operator with entries Vn = 
hn/hn-i, we find the expansions 



k=i 



Q2 = 9a' + Y.^-''^"^^ 

k=i 



(3.5) 



where the s^^\ t^^^ and u^^^ are diagonal operators in the pn basis. For simplicity we shall 
from now on go directly to the planar limit n, ^ 00 as before with x = n/N fixed, in 
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which the operators s^'^\ t^'^^ and u^'^^ become functions of x, while a now plays the role of 
a dummy scalar expansion parameter. The last two lines of allow clearly to express 
the s*^'^-' and t*-'^-' in terms of the u^'^\ Writing moreover the relation Q2 — Q\ — gziQ^ at 
orders 8, 5, 2 in cr we may express u^'^\ u^^\ u^'^^ as 



u 



1 + 2g'^ziZ2v'^ 
= -2g^z^Z2U^'^v' ^3.6) 

We finally get two equations determining u^^^ and v implicitly in terms of x, by writing 
Qi = Q3 + gz2Q\ at order 1 in cr and Pi/N = Q2 — gQ\ at order —1 in a. Upon defining 
the rescaled quantities 

a = — 

Z2 

V = g^Z2V (3-7) 

U = Z2{u^'^f 



we end up with the two equations 

g'zlx = ^(y, U) ^ W^l - 2a'V') - ^V'^^^^^ + ^(1 - 20a'V')U 

. , {l-2aV^ -AUa'^V^) 



Z2 = i;{V, U) = U [ 2V{1 - Aa^V^) + 



[l + 2aV^) 



i3A 



This system generalizes ( ^.29 ) in the sense that we must first solve the second equation 



for U{V) as an implicit function of V (namely Z2 = i{j{V,U{V))) and plug it back into 
the first equation to get the relation g'^z'^x = f{V) = f{V, U{V)), leading to the formula 
for the planar free energy through relation (|2.15|) (upon the substitution g ^ g'^Z2 in the 
denominator of the Log) . More precisely, the correct determination of U is dictated by the 
small X limit in which f ~ a;, hence V ~ g'^Z2X and U ^ Z2- 

Before we turn to the general study of the critical lines of the model, it is instructive 
to analyze the simple limiting cases discussed in Sect. 3.1, namely ^ (^i — > 0) for 
which we expect to recover the one-particle model of Sect. 2.3, and (I/22) (^i, Z2 00 
with u = Z2/z\ fixed) for which we expect a first order transition. 

For Zl 0, we simply take a = in ( |3.8| ) to write the second equation as [/ = 
^2/(1 + 2VY while the first equation gives g'^z'^x = VU — 21/^(1 — 21/)^. We therefore 
recover eqn. ( |2.29| ) with Z2 ^ z. 
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More interestingly, in the other hmit, we must let the parameters scale as zi = zi/e, 
Z2 = -22/^^, g = with e ^ 0. From ( |3.7| ) we deduce the other rescalings a = ae, V = V 
and U = U/e^. In this limit, eqn. (|3.8|) becomes 

g'zlx = l{2UVa^){l - {2UVa')) 
Z2 = U{2V + 1 - 4UV^aY 
The second line of (|3.9| ) may be recast as 

= -{2UVa^) (1 - {2UVa^)) (3.10) 
with = zf/z2 = zf/z2 as before. This gives an alternative expression for 0{V) = g^zfx. 




The maxima of 0{V) correspond clearly to either 2tJVa'^ = 1/2 or \JU /z2 = 1/2 leading 
respectively to the critical values g1 = l/(8if) or = /{9>zl) = 1/(852)- As before, 
the choice of the correct determination is best seen by considering the critical lines in the 
plane {y,u~^). Using ( p.lO| ) with 2tJVa'^ = 1/2, we obtain the first curve 

u = — ■j^=^= (3.11) 

2^Z2Va^ - 1 



Using (|3.10|) with now \jlJ Iz2 = 1/2, we get the second curve 

u"! = z2Vd?{2 - Z2V6?) (3.12) 



The curves ( |3.11|) and ( |3.12|) are plotted in Fig. 8. The solid portions correspond to the 
smallest values of V for fixed which define the location of the relevant maxima of (piy^ 
attained by the change of variables x -^V . The transition between the two curves takes 
place at w = 1 where cp' changes expression from g\ to g\ (see Fig. 9). This is clearly a 
first order transition as the slope of g^{u) has a discontinuity at w = 1. Note also that 
(p" {y^ is non-zero at the transition point, therefore we have 7 = —1/2 as in the case of 
pure gravity. 

Let us now turn to the general analysis of the complete phase diagram as obtained from 
eqn. ( p.8| ). As before, we first look for critical lines characterized by iy^ = = dy^p — 
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Fig. 8: Critical lines in the (F = Z2V6l^,u ^) plane corresponding to 
(f'{V) = 0. The correct line corresponds to the lowest value of V and is 
represented by a solid line. The first order transition point w"^ = 1 is char- 
acterized by the contact between the two determinations of at which 
jumps from 1/(852) to l/(8if). 

du<p{dvip / dutp) ■ Using the explicit expressions for (p{V, U) and V'(^) U)^ we immediately 
get two possible conditions: 



= (1 - 2aV'^f{l - 2V{1 + 2aV'^f) + AUa'^V^{l + 6aV^) 
= W(l + 2V(1 + 2aV^ff 

- U(l - 8aV(l -V- 15aV^ - 52a^V^ - 50a^F^)) + AU^a^V 



(3.13) 



which leads to the three determinations: U = Uq solution of the first line; U = U± 
conjugate solutions of the second line. 

When substituted into Z2 = ipCV^U), this gives three branches in the plane (V, -22) 
for fixed a. Let us now restrict ourselves to -21,-22 > 0, hence a > 0. For small enough 
positive values of a the branches in the {V, Z2) plane take the generic form displayed in 
Fig. 10. As before, the correct solution corresponds to the smallest value of V for fixed 
^2, represented by solid lines in Fig. 10. We clearly identify two transitions at values Z2'^^ 

(c) 

and Z2 ■ Comparing the qualitative behavior of the curves in the vicinity of the transition 
points with the behaviors obtained so far in the two above limiting cases, we can identify 
z^^"^ with a first order discontinuous transition point and Z2^^ with a continuous (critical 
Ising) transition point. 
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1 u 

Fig. 9: Critical line in the {u~^, g'^zf) plane for the limiting case ^1,22 —>■ 00, 
= z\l Z2 fixed. The behaviors of <^(t^) below, at and above the transition 
are displayed in the small plots. The slope of g^z\ is discontinuous at the 
first order transition point = 1. 

The transition points correspond as illustrated in Fig. 10 to the coincidence of two 
branches. We can therefore obtain their location by expressing that = C/+ or C/q = t^-, 
which is done explicitly by writing that the solution t/o of the first line of ( p.l3|) also 
satisfies the second line. We get the two possible conditions: 

= 1 - y(l - 2V'^a - 20(y2a)2 - 2My'^af) 

0=1- 2W'^a + 100(^2^)2 _ 2V{\ - 16V^a - 104(^2^^)2 - 448(^2^)3 + 400{V^a)^) 

(3.14) 

Introducing the variable W = 2V^a, we can express all the relevant quantities as rational 
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V 



Fig. 10: Critical lines in the {V^ Z2) plane for 2:2 > at some fixed (small 
enough) value of a = Zi/z2., as obtained by setting ^'(V) — 0. The solid 
lines correspond to the correct branches given by the smallest values of V . 
They correspond successively to the solutions Uq and Uj^ oi eqn. ( p.l3|) . 

A continuous transition point occurs at the crossing point Z2 = of the 
two lowest branches, with a situation analogous to that found in Fig. 4 at the 
point z = Zs^. A discontinuous transition point occurs at the contact point 



Z2 



between the two upper branches, with a situation analogous to that 



of Fig. 8 at the point u =1. 



fractions of W. For the first line of ( |3.14|) we end up with 

1 W^{1 + W)^{1-3W)^ 



Z2 
1 

u 



4(1 - Wf 
{1-W)^{1 + W)'^ 



(3.15) 



1-3W 

For the second line of ( p.l4|) we get 

1 {1-8W- 26W^ - 56W^ + 25Vr4)5 



Z2 32(l-M^)6(l-5W)4(l-8Vr-25Vr2)2 

1 _ 1281^2,^1 - Wf{l -8W - IhW'^f 
u 



(3.16) 



(1 - 5W)^{1 -8W- 26W^ - 56W^ + 25W^) 
It is easy to check that the points of the second curve ( p.l6|) correspond to crossings of 

(c) 

branches such as that happening at Z2 in Fig. 10. The situation around these points is 
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totally analogous to that described on Fig. 3 at the Ising transition point, with in particular 
f"{V) = 0, henceforth 7 = —1/3. We identify this curve with a line of continuous critical 
Ising transition points (CFT with c = 1/2 coupled to 2D quantum gravity). 

The first curve (|3.15| ) on the other hand corresponds to a contact of branches such as 
that encountered at Z2'^\ The situation around this point is now similar to that found in 
Fig. 9 at the first order transition point. On this line we have f"{V) 7^ 0, hence 7 = — 1/2 
as in the case of pure gravity (CFT with c = coupled to 2D quantum gravity). However, 
the critical parameter (7^ has a discontinuity in its slope across this line, hence so does the 
free energy. We thus identify this curve with a line of first order transition points. 



1st order 




Fig. 11: The phase diagram of the two-particle exclusion model on ran- 
dom vertex bicolorable planar lattices in the {l/z2,u~'^) plane for Z2 > 0, as 
obtained from the exact solution of the six-matrix model ( |3.1| ). This phase di- 
agram agrees with that of Fig. 6 with an Ising-like critical (solid) line meeting 
a first order (dashed) one at a tricritical Ising point (t). We have indicated the 
relevant branches of U solving eqn. ( |3.13|) , namely Uq in the ordered phase, 
and U± in the fiuid one. We have also indicated a typical constant a = Zi/z2 
hyperbola along which we encounter successively the two continuous and dis- 
continuous transitions of Fig. 10. 

These results are summarized in Fig. 11 where the first order transition line is repre- 
sented by a dashed line and the continuous one by a solid line. These two lines meet at a 



27 



tricritical point with 



W^^ = 

5 

1 15633086927- 2441464587^41 

70= mmm = 0.00107319... (3.17) 



32(-664779 + 87849V4T) ^ ^ ^^^^^ 



u 



(*) 78125 



Note that Fig. 10 corresponds to a typical smaU enough constant a section (hyperbola 
in the {1/ Z2,u~^) plane) such as that represented in dotted line on Fig. 11 which crosses 
the two transition lines successively. The ordered phase lies below the transition lines and 
is described by the solution Uq of the first line of ( |3.13| ) while the disordered phase lying 
above the transition curves corresponds to U± solutions of the second line of ( |3.13|) . 

The above phase diagram gives an explicit realization of that described qualitatively in 
Sect. 3.1. To complete our study, let us now show that the tricritical point (|3.17| ) displays 
the expected behavior for a tricritical Ising transition point (CFT with c(4, 5) = 7/10) 
coupled to gravity. A first evidence comes from the fact that the string susceptibility 
exponent at this point is 7 = —1/4 as if' = (p" = (p'" = at this point while cp^'^^ 7^ 0. 
Note that the vanishing of (f" holds generically for all points of the critical Ising line 
( pi.l6| ), as one readily checks by direct calculation. The vanishing of (f'" holds only at 
the tricritical point, as may be checked by a direct calculation too. The tricritical Ising 
CFT is the only unitary CFT with 7 = —1/4 when coupled to gravity!. A second check 
can be performed by computing the thermal exponent a. More precisely, one can define 
two thermal exponents pertaining to two thermal operators with conformal dimensions 
^33 = 1/10 and /i32 = 3/5, and the corresponding dressed dimensions when coupled to 
gravity A33 = 1/4 and A32 = 3/4. As explained in Appendix A, the most relevant one 
($33) governs the generic approach to the critical point through / ~ (22 — ^2*'*)^"" with 
a = {1 — 2A33)/(1 — A33) = 2/3, while the other operator ($32) governs the fine-tuned 
approach along a line tangent to the critical curves ( |3.15| ) and (|3.16| ), with a behavior 
/ ~ (Z2 - 4*^)^""' with a' = (1 - 2A32)/(1 - A32) = -2. To obtain the value of the first 
exponent in our model, a simple procedure consists in first fixing the ratio z\l z^-, expanding 
V for both lines of ( p.l3|) in terms oiV — F*-*) and substituting the result into (|3.8|) . We 



^ 7 = —2/{p -\- q — 1) for a central charge c(p, q) hence p + g = 9 and g — p = 1 by unitarity 
yield p = 4, g = 5. 
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finally get - (^W)^ = a{V - V^^^f + b{V - y^*))^ + ... and ^2 - 4*^ = a'{V - y^*))^ + 
_X/(*))4 + _ _ _ ^hich upon inversion leads to g^^-{g^*^)^ = a" {z2- z'i'^)+h" {z2- zff/^ 
with values of a" and h" 7^ independent of the determination of U , hence 2 — a = 4/3as 
expected. To get the second exponent, we may approach the tricritical point by traveling 
along the transition lines ( |3.15| ) and ( |3.16| ). For the first transition line we have 

2 W'^{l + Wf{l-?,Wf{l-2W + 'iW'^ + 2QW^+?,W^-bQW^-?,hW^) 

32(1-W^)i2 ^3.18) 

For the second transition line we have 

2 _ 15(2 - 12W - 41W^2 ^ 54^3 _ [l-m - 2QW^ - 56W^ + 25]^^)^ 

~ 8192(1 - Vr)i2(l - 5W)^{1 -8W- 25W^2)4 (3. 9) 



Using the corresponding parametric values of Z2 ( |3.15|) and (|3.16|) respectively, we easily 
check that the first three derivatives of g"^ with respect to Z2 match at the critical point 
W = W^^\ while the fourth one is different. This discontinuity corresponds to 2 — a' = 4, 
hence a' = —2. 

For completeness let us finally discuss the case zi > and Z2 < {u~^ < 0). Looking 
again for the critical lines with (f' = 0, the relevant branch for U yielding the physical 
determination of is given by the second line of ( |3.13| ) i.e. U = say. Once substituted 
back into the second line of ( |3.8| ), this gives a critical line in the {V, Z2) plane, such as that 
plotted in Fig. 12 (a) (for a = Zi/z2 small enough). We recover the Lee- Yang singularity 
point at some ;2^(w), characterized by ^p"{V) = (7 = — 1/3) due to the merging of the 
maximum of f with its further minimum. Eliminating U, we finally obtain the following 
parametric curve 



1 _ (5 + 100M^ + 326M^2^ 820M^3_ 2275M^4)^ 
7^ ~ ~ 512(1 - 5W)^{1 + 35W)^{1 + 20W + 35Vl^2)2 

1 2048W^2(i _ 5vt/)6(i ^ 2OW + 35W^)^ 

u ~ ~ (1 + 35H^)4(5 + lOOVT + 326M^2 ^ ^2QW^ - 227bW^) 



(3.20) 



with the parameter W = —aV"^ > 0. In Fig. 12 (a), we note the existence of a further 
maximum of ip never attained by the change of variables x ^ V except when these three 
extrema merge simultaneously (see Fig. 12 (b)). This point corresponds to a higher order 
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Fig. 12: Critical line in the (V, Z2) plane for 2:2 < as obtained by setting 
<fi'{V) — for fixed (small enough in modulus) a — Zi/z2. As usual, the 
correct portion of the curve corresponds to the lowest value of \V\ and is 
represented by a solid line. The Lee- Yang critical point Z2 is characterized 
by dz2/dV = and corresponds to the merging and annihilation of the two 
extrema 1 and 2 in (a) with (p"{V) = 0. By increasing \a\, we reach the 
situation (b) where a third extremum 3 merges with them. This defines the 

tricritical point Z2 satisfying in addition ip"'{V) = 0. 

multicritical point with ip'" = 0, hence with a string susceptibility exponent 7 = —1/4. 
This tricritical point corresponds to the values 

455 

1 172647361044v/T4 + 645414154777 , 

^TTT = = -4.07373... (3.21) 

317007031250 ^ ^ 

1 _ 524288(2401452v^ - 8699159) ^ _^ quqqy2 



u 



10274243531825 



We thus conclude that the line of Lee- Yang critical points Z2 (u) ends at this higher critical 
point Z2 ■ We identify this point with the only non-unitary CFT with 7 = —2/{p+q—l) = 
—1/4, i.e. c(2, 7) = —68/7, coupled to 2D quantum gravity. A confirmation of this fact 
may be obtained by computing the thermal exponents of the theory a and a' characterizing 
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the singularity of the free energy (or equivalently of g^) as Z2 approaches the tri critical 

it') 

value Z2 ■ Like in the unitary case, the exponent a governs the generic approach to 
this point in the (1/2:2, plane, while a' governs the fine-tuned approach along the 
line of Lee- Yang singular points. As shown in Appendix A, the predicted values are 
a = 2/3 and a' = 1/2. As before, we obtain the first exponent in our model by expanding 
and Z2 along a curve with, say, 21/22 =const. We find generically — ((7*-* '')^ = 
aiy - y + h{V - y + ... and 22 - 2^ ^ = 0! {V - V^'"^f ^h' {V - V^^'^f + . . . which 
upon inversion leads to g1 — {g'^^'^Y = '^"(^2 — ^2* ^) + ^"(^2 — ^2* ^)^^"^ with h" 7^ 0, which 
gives 2 — a = 4/3. To compute the second exponent, we use the parametric equations 
( p.20|) and the corresponding value of g^ 

2 _ 15(5 + + 326W^^ + 8201^3 - 2275W^^)^(2 - 601^ + 151M^^ + 630W^^ - 2275M^^ 

~ 2097152(1 - 514^)12(1 + 35Vl^)2(l + 2m + 35M^2)4 

(3.22) 

and expand (7^ and 22 around the tricritical point W^^ \ with g\ — [g^^ •')^ = aiy/ —W^^ •')^ + 
h{W - W^^'^f + ... and 22 - z^^^ = a'{W - Vf^(*'))2 + b'{W - Vf^(*'))3 + . . . which upon 
inversion leads to g'^ — {g^^"'Y = a" (22 — -^2* + ^"{^2 — Z2 ''))"^/2 with h" 7^ 0, which gives 
2 - a' = 3/2. 

As discussed at the end of Sect. 2.2, the Lee- Yang critical line separates a small I22I 
region where gc is real from a large I22I region where it becomes complex and generates 
an oscillatory behavior of the form ( |2.2C1| ) for the canonical partition function ^^(21,22) 



for planar graphs of fixed area A. This separation between oscillatory and non-oscillatory 
behaviors extends beyond the tricritical point (t') in the form of a first order transition 
line as depicted on Fig. 13. This line corresponds to a situation where the two complex 
conjugate values of g corresponding to the complex conjugate solutions to = cross in 
modulus the real value of g corresponding to the real solution to (p' = 0. The transition 
through this line is first order in the sense that the thermodynamic free energy has a 
discontinuous slope across the line. The line clearly originates at the tricritical point {t') 
where the three values of g are real and coincide. More interestingly, it is easy to see 
that it terminates at the point {l/z2,u~^) = (0,-1). Indeed, along the axis 1/^2 = 0, 
we have to compare the real value gi = 1/ \/8£i^ and the two complex conjugate values 



§2 = ±z/ a/8 1 22]. For > —1, we have an oscillatory behavior of the partition function 
Za{zi,Z2) = lime'^Z'^(2i, ^2) ~ (—8^2)^^^, while for < —1 we simply have a non- 
oscillatory behavior Za{zi,Z2) ~ (2v^2i)^. 
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Fig. 13: The phase diagram of the two-particle exclusion model on random 
vertex bicolorable planar lattices in the {1/z2tU~^) plane for Z2 < 0, as ob- 
tained from the exact solution of the six- matrix model (3J.). A (solid) line 



of Lee- Yang type critical points terminates at a tricritical point {t') where it 
meets another first order (dashed) transition line. These two lines form the 
border of the region of the phase diagram where an oscillatory behavior is 
observed in the canonical partition function Za{zi, z^)- 

4. Conclusion and Discussion 

In this paper, we have shown how to reach critical and multicritical points in the con- 
text of nearest neighbor exclusion models. We have displayed a number of exact solutions 
for various models defined on random planar lattices. The crucial outcome of this analysis 
is the importance of the colorability condition of the underlying lattice. Indeed, we have 
shown that critical points based on a ^2 symmetry such as the critical Ising and tricritical 
Ising models are reproduced for exclusion models under the condition that the lattice itself 
be vertex-bicolorable. In view of understanding the physics of the corresponding models 
on regular lattices, the outcome of the random lattice solutions is twofold: first it sorts 
out which features of the lattice itself take a relevant part in the models' critical behav- 
ior; secondly it allows for an exact solution in the planar limit that clearly identifies the 
critical universality classes for regular lattice models (such as the celebrated Hard Square 
model) which are still unsolved directly. More precisely, this study allows to infer that 

32 



the bipartite nature of the square and honeycomb lattices indeed translates into a critical 
Ising-like universality class for their crystallization transition point. 

We have shown how to extend the definition of hard particle models to reach higher 
order critical points, by introducing a two-particle exclusion model in which each site may 
be empty, singly or doubly occupied with an exclusion constraint that a total of at most 
two particles may share the same edge. A straightforward generalization consists in having 
a /c-particle exclusion model, where sites can be occupied by 0, 1, 2, ... up to k particles, 
with a weight Zi for an occupancy by i particles, and with the exclusion constraint that a 
total of at most k particles may share the same edge. When defined on a regular lattice, 
the exclusion constraint is easily turned into an non-overlapping constraint for appropriate 
tiles of various sizes. When defined on a random vertex-bicolorable lattice, these models 
are amenable to a {2k 2) matrix integral, with potential 

2k+l k 
V{A,, A2, A2k+2) = Yl (-1)^+'^.^.+! - I E ^^(^2.+l + ^2(fc-0+2) (4.1) 

for trivalent graphs, with a weight g per vertex (and zq = 1). Remarkably enough, the 
exclusion rule translates into a chain-like quadratic interaction between the matrices, which 
makes the models exactly solvable by means of standard orthogonal polynomial techniques. 
With the model ( [4.1j ), we expect to be able to reach multicritical points governed by CFTs 
with central charges c(p, q) with p + q = 2k + 3, hj appropriately fine-tuning the z^'s, in 
order to reach the highest possible critical string susceptibility exponent 7 = — + 2). 
We also expect the situation to be identical on vertex-bicolorable regular lattices. 

One may also address the question of the behavior of the same particle exclusion 
models on lattices with other colorability properties. From the exact solution of the Hard 
Hexagon model on the triangular (hence vertex-tricolorable) lattice, it is natural to expect 
that the above exclusion models, when defined on vertex-tricolorable (fixed or random) 
lattices, give rise to critical and multicritical three-state Potts models. Unfortunately, 
the corresponding matrix models are no longer directly solvable by means of orthogonal 
polynomials. 

More generally, the vertex-/c-colorability of the lattice is likely to generate within the 
context of exclusion models critical points with symmetry such as the /c-state Potts 



model or the 2^ model of Fateev and Zamolodchikov [^. We know however that for 



large enough k > 4 we lose the continuous transition in the Potts case, while in the other 
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case the central charge c = 2{k — l)/{k + 2) > 1 forbids any meaningful coupling to 2D 
quantum gravity. 

Note finally that the above study provides an example of critical phenomena whose 
coupling to gravity is sensitive to the type of graphs summed over. While the precise 
connectivity (tri- or tetra- valency) of the lattices is unimportant, their Eulerian (bipartite) 
or Euclidean (arbitrary) character is relevant. This situation is reminiscent of that of fully- 
packed loop models pi| , P^ . 
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Appendix A. String susceptibility and thermal exponents for minimal CFT 
coupled to 2D quantum gravity 

The matrix models in general give a discretized version of matter systems coupled to 
2D quantum gravity, in the form of discrete statistical models defined on random graphs 
accounting for the fluctuations of space-time. In the continuum approach to 2D quantum 
gravity, one may relate the properties of some critical matter system in fixed geometry 
(typically a CFT with central charge c in the plane) to that of the same matter system 
coupled to 2D quantum gravity. The precise connection involves the celebrated KPZ 
formula expressing the string susceptibility exponent for a unitary CFT of central 
charge c coupled to 2D quantum gravity: 

7(c) = (A,) 

as well as that the dimension A of the gravitational dressing (p of a primary operator of 
the CFT with dimension h: 



A(h.c)^ ^^^Il!^-^ (A.2) 

This dimension measures the singular behavior of the corresponding two-point correlator 
as {(jxp) ~ — iJ.)'^^^'^'^'^~''^'^\ when the cosmological constant n approaches its critical 
value He- The above formulae simplify drastically when considering minimal CFT with 
central charge 

c(p, q) = l- eil^ (A.3) 
pq 
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with say q > p + 1 and p A q = 1, and operator spectrum 

^ = iSI-ZP^l^SP^ (A.4) 
Apq 

with qr — ps > 0, l<r<p — 1, l<s<g — 1. Eqns. ( |A.1|) and (|A.2|) reduce in this case 
to 

7(c(p, g)) = 1 - ^ A,,, = A(/..,„ cfe q)) = ir-l)q-is-l)p ^^ ^^ 

p 2p 

For a unitary theory {q = p + 1), the most relevant operator of the theory is the 
identity with /ii i = and its gravitational dressing known as the "puncture operator" 
P has Ai^i = as well. This operator measures the area of the random surface, which 
allows to identify the deviation Hc — fJ- with Qc — g in the corresponding matrix model, 
and to interpret 7 = —1/p as the exponent governing the singularity of the free energy. 
In the case of pure gravity (p = 2, g = 3), we have 7 = —1/2, and the only operator is 
P. In the case of the critical Ising model (p = 3, g = 4) we have 7 = —1/3 and three 
operators: the puncture P, the dressed energy operator $2,1 and the dressed spin operator 
$2,2, with respectively Ai^i = 0, A2,i = 2/3 and A2,2 = 1/6. The dressed energy is the 
thermal operator, with coupling z^ — z ~ (/Uc — /u)^"^^'^ = (fi'c^fi')^^^? hence upon inversion 
we get a singularity of the free energy of the form given by (^c — 5* ~ (-^c — -2)^"", with 
a = —1. In the case of the tricritical Ising model {p = 4, q = 5), we have 7 = —1/4 and 
there are six operators, among which we distinguish the puncture P and the two thermal 
operators $3,3 and $3,2, with respectively A3^3 = 1/4 and A3^2 = 3/4. The generic 
thermal perturbation of the model is governed by the most relevant operator $3^3 with 
coupling Z2 — Z2^ ~ {fXc — = (g^ — fi*)^^^, yielding the singularity of the free energy 

gc — g ^ (-22 — ^2*'')^"" with a = 2/3. By fine-tuning the parameters one may approach the 
tricritical point on a line along which the contribution of $3^3 is cancelled, and therefore 
the next most relevant thermal operator $3 2 takes over. This leads analogously to another 
thermal exponent a' with 2 — a' = 1/(1 — A3^2), hence a' = —2. 

For non- unitary theories with q > p + 1, the above formulae must be interpreted 
carefully to account for the fact that the identity is no longer the most relevant operator. 
Indeed, the scale of the deviation from criticality is set instead by the operator of smallest 
(negative) dimension ho = (1 — (p — q)^)/{4:pq) (with qr — ps = 1) corresponding to 
the gravitationally dressed operator $0 with dimension Aq = {p ~ q + l)/{2p). The 
latter operator has the coupling — which allows to identify the deviation from 

criticality {gc — g) in the matrix model as gc — g = {fJ-c — A^)^~^°- Note that this general 
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relation also holds in the unitary case q = p + 1, where it reduces to gc — g = fic — IJ', as 
ho = Aq = 0. The most singular part of the free energy is due to the presence of $o and 
may be obtained by writing that (i'^//(ifi'^|smg ~ (^o^o)- We can therefore write 

($o$o) ~ (^c - ^)2Ao-7(c(p,,)) ^ _ gy-r (A.6) 
from which we deduce the corrected string susceptibility exponent 7 of the matrix model: 

' 1-Ao P + q-1 ^ ' 

To compute the thermal exponent we need to identify the next most relevant thermal 
operator say $1 with dimension Ai with coupling proportional to — ~ {gc — 

^^i/(2-a)^ leading to 2 - a = (1 - Ao)/(l - Ai). 

For the Lee- Yang edge singularity, we have p = 2, q = 5, and therefore Aq = Ai 2 = 
— 1/2 and 7 = —1/3. In this case, the deviation from the critical "temperature" Zc — z 
is coupled to the next most relevant operator of the theory, which turns out to be the 
puncture operator $1 = P with Ai = 0, leading to a = 1 + Aq = 1/2. For the case 

p = 2,q = 7 of Sect. 3.2, we have Aq = Ai^s = —1 and 7 = —1/4. A generic deviation 

(t') 

from criticality Z2 — -^2 is coupled to the operator $1 = $12 with Ai = Ai 2 = —1/2, 
hence yields a thermal exponent a = 2/3. In a fine-tuned approach to the critical point, 
we may cancel the contribution of $o? in which case $1 now plays the role of the most 
relevant operator, with Aq = Ai, while the thermal operator becomes P, with A'l = 0. 
We deduce the fine-tuned thermal exponent a' = 1 + Aq = 1-l-Ai = 1/2 identical to that 
of the Lee- Yang case. 



Appendix B. Hard particles on a random tetravalent lattice: double scaling 
limit 

In the following, we derive the double scaling limit of the model of hard particles on a 
random lattice (non-necessarily vertex- bicolorable) . We show that the renormalized string 
susceptibility obeys the Lee- Yang differential equation. We first rewrite the eqns ( |2.11J ) 
and ( |2.12| ) in components, namely 

^ / / N 

— = Vn - Tn - gz[Sn + S^+l + Sn+2 + TniVn-l + + r„+i) 

fn=Vn{l+ gz{rn-l+rn+rn+l)) (B.l) 
Sn = gZVnVn-lVn-2 
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and 

= ^^n - 9{Sn + Sn+1 + Sn+2 + ^'nl^'n-l + ^^n + ^^n+l) 

rn = gvn{fn-i + r„ + fn+l) (B.2) 
Introducing the new coefficients 

Rn = grn, Sn = q'^ Sn-, Rn = gfn, Sn = g'^Sn, K = gVn (B.3) 

we simply get 



Rn 


= Vn{l + z{Rn-l 


-\- Rn + Rn+l)) 


Sn 


= zVnVn-lVn-2 




Rn 


= VniRn-l + Rn 


+ -Rn+l) 


Sn 


= VnVn-lVn-2 




n 


RnRn 

-Vn 

' n 


zVn{Vn-lVn-2 + Vn 



(B.4) 



where both first lines of ( p.l| ) and (p.2|) turn out to be equivalent to the last line of (|B.4|) . 
When N becomes large, all sequences tend to smooth functions of x = n/N , and setting 
a = 1/N, we now make the following scaling ansatz on Vn = V{x) = V{1 — a'^v{x)), 
for some unknown function v{x) for which we will derive a differential equation. We 
must first solve the first and third lines of ( [B.4| ) for r(x) order by order in a, where 
= R{x) = R{1 — a'^r{x)), and where the values of V,R,z are taken along the critical 
line (|2.18|) , namely with 

{l-3u)^ ' 6{l + 3u)' 12(l + 3ti)2 ^ 

with the result 

- 12(1 + 42u + 81u^)v{x)v"{x) - 288w(l + 3u)v'{xf + 432w(l + 3u)v{xf) + O(a^) 

(B.6) 

We now expand the last line of ( [B.4|) up to order 6 in a, after setting u to its critical value 



u 



c 



{2^/5 — 5)/15, with the final result 

^^Llll = a6(^;(a;)3 - -v'(x)^ - v(x)v"(x) + —v^^\x)) (B.7) 
gc 2 10 
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Upon introducing the renormalized cosmological constant y = {Qc — 9x) / (a^^'^gc) and 
appropriately rescaling v{x) a~^^/^f (y), we finally get the standard differential equation 
for the renormalized string susceptibility v{y) 

y = v{yf - \v\yf - v{y)v"{y) + ^v^^\y) (B.8) 

which is easily identified with that of the Lee- Yang edge singularity coupled to 2D quantum 
gravity [@]. 



Appendix C. Hard particles on arbitrary random trivalent lattices 

In the following we show that the trivalent lattice version of (|2.1|) leads to the same 
qualitative physics, namely a unique critical point in the universality class of the Lee- Yang 
edge singularity coupled to 2D quantum gravity. 

We start with the trivalent version of the matrix model (12.11): 



1 53 A3 (^-^^ 

y(A B) = --A'+AB-g—- gz— 

where A, B are Hermitian with size N x N, and the measure is normalized so that 
Z]^ (0, 0) = 1. Comparing with ( |2.1| ), the bi-orthogonal polynomials are no longer 
even/odd as V is no longer even. We may still write the trivalent version of (p.7|) : 

-i = - Qi - gzQj 

p (C.2) 

-^=Q\~9Ql 

and show that the Q's have finite range, with 

Qi = a + r + a~'^s + a~'^t 

(C.3) 

Q2 = r + r + r h + T h 

where a, r denote the shift operators acting respectively on the left and right bi-orthogonal 
polynomials. In components, the equations ( |C.2|) now read 

n , , , 

— = Vn - Sn - gz{tn + tn+1 + Snl^^n-l + Tn) 
rn = {rn+ gz{Sn + Sn+l+rl)) 
Sn = Vn{l+ gz{rn + r^-l)) 



tn = gzvnu 



n-1 



38 



and 

77, ~ ~ 

= 'f^n - gitn + in+1 + Snif^ + f„-l)) 

Tn = 9{Sn + Sn+1 + rl) 
tn = gVnVn-l 

For large n,N with x = n/N, we get algebraic equations for V, R, S,T[x) respectively 
limits of g'^zVn, gzr^i, g'^zsn, g^ztn and their tilded counterparts: 



g'^zx = 


V -S-2zT- 2RS 


R = 


R{l + R) + 2zS 


S = 


V{1 + 2R) 


f = 




R = 


2S + R^/z 


S = 


2VR/Z 


T = 


V^/z 



(C.6) 



easily solved in the form of an algebraic equation for as a function of V 

f R{l + R? \_ 2zV 

and a master equation for the dependence on x: 

9 9 . ^ . ^ V Ril + R) il + 2R) , , 

g\^x ^ ^iy) = zV[\ - 2V) - 2 ^ ^Jl^ ^ (C.8) 

Writing that = = 0, and eliminating z and V ^ we are left with a sixth order equation 
for the critical value of i2, namely 

13 + 266i? + 1810i?2 + 5920i?^ + 10200i?^ + 8748i?^ + 2916/?*^ = (C.9) 

only one root of which leads to a positive value of g'^z^ through (|C.8| ) at x = 1. This leads 
to a unique critical point at 

R^ = -0.090430 ... K = -0.036173 . . . 

(C.IO) 



z^ = -0.16565 ... gc = v/^(K)/kc| = 0.28105 . . . 

with the critical exponent ^str = —1/3, and the critical point is in the (non-unitary) class 
of the Yang-Lee edge singularity as in the tetravalent case. It is easy to check that the 
scaling ansatz = V{1 — a^f (x)), a = l/N, still leads, upon solving (|C.4|) ( p^ ) order by 
order in a, to the same differential equation (|B.7|) in which g g^z'^. 
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Appendix D. Hard particles on trivalent bicolorable graphs: double scaling 
limit 

In the following, we complete the identification of both tricritical points ( p.33|) for 
the model (|2.21|) by deriving the corresponding differential equations for the renormalized 
version of the string susceptibility V{x). To derive the double-scaling limit of the matrix 
model ( |2.21|) , let us first write in components the complete equations (|2.28|) : 

s'^' = 9 

S^n^ = -g^ZVnVn-lVn-2Vn-3 
Vn = + g^ZV^iSn-l + (D.l) 



r. 



_ .(1) _ of^(l) , ^(1) ) 



Setting Vn = g'^zVn, Sn = g'^zsn \ we finally get 

= S'n — Vn{^ — Sn-1 — Sn+l) 

g^z^- = zSnil - Sn-1 - Sn+i) - + Vn+l) (D.2) 

+ Vn-lVnVn+l{Vn+2 + K-2) + K(K-lK-2K-3 + 

We make the following scaling ansatz V{x) = V{1 — a^v{x)), a = 1/N a small parameter. 
We first solve the equation Vn{l — S^-i — •S'^+i) = Sn order by order in a, with the result 
S{x) = S{1 — a^s(x)) 

'^"^^"1 + 2^ (1 + 2^)2'' (1 + 2^)2''^'^^ 12 {1 + 2VY ^ (D.3) 

4 , V{l-AV) 4 

^ (1 + 2^)3 '^ " (1 + 2^)2 '^^^^'^ (1 + 2^)3 '^^'^^ ^ ^ 

which we then substitute into the second line of ( p.2| ) and Taylor-expand in a. Apart from 
the term of order zero, ^{V) = g'^t^ , the lowest order terms are in . We then simply get 
at leading order 

a^{A,v{xf + Biv{x)v"{x) + C,v'{xf + D,v^^\x)) + O(a^) 

A+ = — , A_ = - 
+ 15 3 

= B_ = -2 (D-4) 

+ 5' 

C+ = C_ = -1 



gin 



5 

5' 10 



40 



where the index i = ± refers to the critical point z±. Upon setting y = {gftf — 
g'^t'^x)/{gftfa^/'^) and rescahng respectively t'(x) — > 9v{y)/{4a^'^^'^) andt'(x) 3v{y)/{2a^'^^'^), 
the differential equations take the standard form 

(1) - vv" - + ^v^'^ = y 

(2) v'-vv"-^ivr + {^v^'^=y 

which we immediately identify with the differential equations governing the double scaling 
limit of respectively the Ising model and the Lee- Yang edge singularity M. 



Appendix E. Hard particles on tetravalent bicolorable random graphs 

The study is quite analogous to that of Sect. 2.3. The matrix integral takes the same 
form as ( p.21| ), but with the potential 

V{A,, A2, As, A4) = A,A2 - A2A, + AsA^ - ^{Aj + Aj) -^{At + Aj) (E.l) 

We still have the symmetry Ai ^ A^-i for i = 1,2, hence introducing again a family 
of monic orthogonal polynomials Pn, and keeping notations as n ( |2.23|) , we get the two 
equations 

— =Q2- 9Q\ 

N (E.2) 

Q = Qi-Q\- gzQl 



The potential now satisfies the following additional symmetry property, replacing ( p.24| ): 



— ZX2, ix3, — ZX4) = X2, Xs, X4), = — 1 (E.3) 



The operators Qi,Q2 still have finite range, as a consequence of ( [li).2| ), and read, thanks 
to the symmetry ( [lii.3| ): 

7 



-4jV(i) 

(E.4) 



with the usual shift operator a, with a'^ = a ^f, where Vn = hn/hn-i as usual, hn = 
iPmPn) the square norm of Pn- The adjoint of Q2 reads 

q\ = g{a-\f + s^^^v-^a + s^^\v-^af + s^^\v-^af (E.5) 
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and the equations ( [E.2|) boil down to 

,(3) = -g^za\a-\f 

= - gz{as^^^a-^ + (7-^.(2)^2 ^ ^-5,(2)^5) 

+ + aV(i)a-5s(2) ^ ^-2^(2)^2^(1) ^ ^2^(1)^-2^(2) 

" ,(l)_^(^(l)+^-V(l)^ + ^-2^(1)^2) 



or equivalently in components: 



r. 



'n ^ = -£?^2;t^n^^n-l---t'n-4(si-5 + ^-2 + «n+l) 

^ I O^Co2^ (3) , (3) , „(3) X , . (1) (2) , ,(2) (1) ,.(1)5(2) 



+ .(2)Ji) , Ji) J2) , J2) (1) ^ , ,(1)J1) Jl) ^ 
!}_ _ „(1) _ of^(l) I ^(1) I ^(1) ^ 



(E.7) 



Upon the redefinitions 

V^ = gvr,, R^ = gM^\ = gs^^\ Tr, = gh(^\ Ur, = gh^^^ (E.8) 

these equations finaUy reduce to 

Un = -zVnVn-l...Vn-8 

Tn = —zVnVn-l---Vn-A{Sn-b + Sn-2 + 5'n+l) 

'S'n = Vn{l — 2(Tn,-i + + ^n+S + SnS^-l + S'„-i5'^_|_i + S'n+l5'„+2)) 

= VnVn-lVn-2 + ^(?7n + t^n+3 + 

+ SnTn-1 + TnSn-5 + SnTn+2 + TnSn-2 + 5'n+3^n+2 + 7n+3'S'n-2 + 5'n5'ri-l5'n-2) 
= 'S'n — (-Rn + -Rn+1 + -Rn+2) 

(E.9) 
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In the planar limit, each sequence tends to a function of x = n/N, n,N ^ oo, which 
we label by the same capital letter. Namely writing U = —zV^, T = —3zV^S, and also 
introducing E = S/V we find 

1 = E{1 - 9z^V^) + 3zV^T.'^ 
gx = ip{V) =EV- 3{V^ + zV^E^ - ISz^V^E"^ - 3z^V^) (E.IO) 
= -31/3(1 - 3z'^V^) + 1/E2(1 + 45^2^6) 

where we have used the first equation to simplify the second. The critical line is the 
solution of ip'{V) = 0. The first line of ( [bil.lO| ) allows to compute 

_ 6zVE{9zV'^ - S) 
dV ~ 1 + QzV^T, - 9221/6 ^ ■ > 

Setting W = zV^ , we find 

_ ((1 - bW)T. - W{1 + W)){{1 + 5Vr)E + 3^(1 - W)){2WT. - V{1 - W^)) 
~ {2WT, + V{l-W^)) 

(E.12) 

hence we have three solutions 

^ _3V(1-W) 

^ ^ 1+5W 

(3) E=^'^-'^'' 



2W 

The first line of (FTOD allows to express V in terms of W: 

1 fl-5W 



(1) V 



3(1 - 3M/ + 5M/2) \1 + W 



(2) V = f l±^f (E.14) 

^ ^ 3{1 + 3W + 5W^) \1-W ' 

(3) V 



We deduce the value of 2; = 



(1) z 



3(1- Vr2)2 

91^(1 + W)^{1 -3W + 5Vr2)3 
(1 -5W)6 



^ 9W{l-Wni + 3W + bW'r 
^ ' (l + 5W/)6 

_ 9(1 
^~ 64W^ 
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and finally that of g = ^p{V): 

2(1 - 5W)^{3 + 24W - 101^2 ^ 4Qy^3 ^ 35^4^ 



:i) 9 



27{1 + W)^{1-3W + 5W^)^ 



_ 2(1 + 5Wy{3 - 24W + lOW^ - 40W^ + 35W^) 

^ ' ^~ 27{1-W)^{1 + 3W + 5W^)^ ^ ■ ' 

16W(1 + W^)(l - WW + 5Vr2) 
(3) g = 



27(l-M/2)6 

The tricritical points are solution in addition of f"{V) = 0. In the three above case, we 
get 

{18V{l-12W+5W'^){-l-20W+35W'^) ■ /.n 
(T35vi7)2 m case [L) 
lSVil + 12W+5W^)i-l+20W+35W-) ^^^^ ^3) (E.17) 
-^1^(1 + 5W^){1 - 12W + 5W'^){1 + 12W + dW^) in case (3) 

As in the trivalent case of Sect. 2.3, each of the lines (1) and (2) have two tricritical points, 
one of which is a cusp, the other coming from the intersection with the critical curve (3). 
Moreover the curve (1) is always reached before (2). The cusp solves 1 + 20W — 35W^ = 0, 
while the intersection solves 1 — 12W + 5W^ = 0. Again, one of the two branches of 
these equations is always reached before the other. We end up with a qualitative picture 
identical to that of Fig. 3, with the cusp and intersection corresponding respectively to the 
values W±: 



10-3V15 6-V31 , 
W- = W+ = ^ (E.18) 

while the tricritical point corresponds to the values {z±,g±) given by the case (1) in ( [E.16| ) 
and ( [E.15| ), respectively with W = W± of ([E.18|) , with the exact values 



15683(-83151 + 26080v/l5) ^^^^^^ 

= — .13d5d8... 



2573571875 



2000 129105 - 24881 V15 „„„„„ , 

o_ = ^ — ^ = .277724... 

^ 235782657 



^ 6561(146327- 21472^31) ^ 

+ 9765625 



(E.19) 



^ 80(14903 - 1067V/31) ^ 
14348907 

To get the double scaling limit, we set 

Vn = V{l-a^v{x)), Sn = ^V{l~a^s{x)) (E.20) 
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and first solve the third equation of (|E.9|) order by order in a, with the result 



4 Vr(l - 5M^)(3 - 155W^ + 650W^ - 1530W^^ + 13751^^ - 1751^^) 3 
^ W{1 - 5W){3 - 205W + 355W^ - 85W^) 2 

+ (1 + W^)3(l + 5W^2) ^ 

4l2M^(l-5M^)(l-29M^ + 70M^2_i90^3^225M^4-25Vr5) , , , 
+ " (1 + H^)3(1 + 5M^^)^ 

(E.21) 

We then substitute this and (|E.20| ) into the last line of (|E.9| ), and Taylor-expand in a. The 
final result reads 



9i 



- gx 

h 




A_ 






2 


B_ 


= -5D. 


C- 






2 




27(5 - 



a^{A,v{xf + B,v{x)v'\x) + dv'^xf + DiV^^\x)) + O(a^) 

A+ = -D+ 
+ 2 + 

B+ = -\d+ (E.22) 
9 

= -Id, 

r. _ 291 - 49^31 



110 ^ 5125 

which may be put back in the standard forms ( p.5| ). We conclude that the situation 
for the hard particle model on vertex-bicolorable tetravalent random lattices is qualita- 
tively the same as that for trivalent ones: we find two critical points at Z- = —.136568... 
and Zj^ = 17.989334... respectively in the universality classes of the Lee- Yang and Ising 
critical models on random lattices. Note that this model is precisely the gravitational 
version of the classical Hard Square model, as, in the dual picture we are considering 
non-overlapping square tiles on random, face-bicolorable planar quadrangulations (dual 
to planar tetravalent vertex-bicolorable random lattices), hence a random version of the 
square lattice incorporating its vertex-bicolorability. 
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